A Molecular Matter- Wave Amplifier 
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We describe a matter-wave amplifier for vibrational ground state molecules, which uses a Feshbach 
resonance to first form quasi-bound molecules starting from an atomic Bose-Einstein condensate. 
The quasi-bound molecules are then driven into their stable vibrational ground state via a two- 
photon Raman transition inside an optical cavity. The transition from the quasi-bound state to the 
electronically excited state is driven by a classical field. Amplification of ground state molecules is 
then achieved by using a strongly damped cavity mode for the transition from the electronically 
excited molecules to the molecular ground state. 

PACS numbers: 03.75.-b,03.75.Pp,42.50.Pq 



Recent progress in the application of Feshbach res- 
onances to ultracold atomic gases has led to the pro- 
duction of macroscopic numbers of ultracold molecular 
dimers starting from either a Bose-Einstein condensate 
(BEC) lj or a degenerate Fermi gas |2|. This work has 
culminated in the production of molecular Bose-Einstein 
condensates starting from a Fermi gas of 40 K [3( and 
6 Li 0. The ability to coherently produce a macroscopic 
number of quantum degenerate molecules opens up new 
avenues of research in the held of matter-wave optics 
such as molecular interferometry, the molecular laser, 
and nonlinear wave mixing between atomic and molecu- 
lar fields. In fact, the effective Hamiltonians that describe 
the coherent atom-molecule conversion are formally the 
same as second harmonic generation in nonlinear optics 
with a nonlinear susceptibility. 

A major stumbling block on the road to coherent 
molecular optics is the short lifetime of the molecules 
formed via Feshbach resonances. While these molecules 
are translationally very cold, they are vibrationally hot 
and can decay to lower lying vibrational states via in- 
clastic collisions with other atoms and molecules. For 
an atomic BEC, the molecules are expected to decay 
in a time ~ 10 — 100/zs for typical atomic densities 
Q. For atomic Fermi gases, recent results indicate that 
the molecular lifetime can be several orders of magni- 
tude longer due to the underlying Pauli blocking between 
atoms 

Still it would be desirable to create molecules in their 
vibrational ground state so as to maximize the amount 
of time with which experiments can be done with them. 
This can in principle be done using two-photon stimu- 
lated Raman photoassociation 0. This is a relatively 
inefficient process due to the very small Frank-Condon 
factors between the free atom pairs and the molecules. 
This is a consequence of the small spatial overlap between 
the tightly bound molecular wave function and the large 
average separation between pairs of atoms in these di- 
lute gases. One way to overcome this difficulty is to use a 
Feshbach resonance to form weakly bound molecules, fol- 
lowed by a two-photon Raman transition to transfer the 



molecules to a low lying vibrational state |8| ■ The advan- 
tages of this procedure are two- fold. First, recent exper- 
iments with Feshbach resonances in quantum degenerate 
gases have demonstrated the ability to readily convert 
a significant fraction of the atoms into molecules. Sec- 
ondly the spatial overlap between the resulting weakly 
bound states and the ground-state molecules is substan- 
tially larger than for that of the free atom pairs, which 
results in larger two-photon transition rates. Indeed, 
several experiments using single-photon photoassociation 
spectroscopy have observed a strong enhancement of the 
photoassociation rate in the presence of a Feshbach res- 
onance |9j. 

In this paper we use the idea of Feshbach stimulated 
photoproduction || as the basis for a matter-wave am- 
plifier for ground-state molecules. Active elements such 
as matter-wave amplifiers are an essential element for 
the field of matter- wave optics and have previously been 
demonstrated experimentally for atoms |lCj . 

In our scheme, pairs of atoms in a BEC are coupled 
to a quasi-bound state of a closed collisional channel via 
a Feshbach resonance • A two-photon transition cou- 
ples the quasi-bound molecular state to the vibrational 
ground state of the molecules. A large amplitude classi- 
cal electromagnetic field drives the the quasi-bound state 
to an electronically excited molecular state, followed by 
the emission of a photon into a single mode of a quantized 
optical cavity field. This second step takes the molecule 
from the excited state to its vibrational and electronic 
ground state. By using a strongly damped cavity one can 
achieve continuous unidirectional coherent amplification 
of an initial molecular signal field that persists until the 
atomic condensate is depleted [l^. Although the num- 
ber of photons in the cavity field is less than one at all 
times, strong coupling of the molecules to the cavity field 
can be achieved by controlling the volume and geometry 
of the cavity as was recently demonstrated in the optical 
regime with the single atom laser 0] . 

The discrete mode structure of the cavity allows one 
to select a particular vibrational state of the electronic 
ground-state molecules provided the cavity linewidth, k, 



2 



is less than the vibrational level spacing, 8v ~ 1GHz. 
This is a distinct advantage over using free space sponta- 
neous emission to initially populate the ground state, as 
in atomic matter-wave amplifiers 10], since spontaneous 
emission would populate a large number of vibrational 
levels of the ground state as well as leading to dissociation 
back into the continuum 01 . As a result, a molecular 
matter-wave amplifier relying on free space spontaneous 
emission would require that the initial signal field that 
is to be amplified be very large in order that stimulated 
scattering dominates over spontaneous decay into other 
bound vibrational and continuum states. 

Atoms in the BEC are assumed to all occupy the same 
hyperfine state denoted by |0). Pairs of atoms in |0) 
are coupled to the quasi-bound molecular state |1) via a 
Feshbach resonance with coupling strength a and with 
an energy difference between pairs of zero energy atoms 
and the quasi-bound state given by the difference in the 
Zeeman interaction energies for the two states, lji 
State |1) is coupled to an electronically excited molecu- 
lar state, 1 2), via a classical optical field with Rabi fre- 
quency rti and frequency u>i. State |2) is also coupled 
via a single mode of the cavity field with vacuum Rabi 
frequency g and frequency lu c to molecules in their elec- 
tronic and vibrational ground state, |3). The internal 
energies of states |2) and |3) relative to pairs of atoms 
in |0) arc lo^ and w 3 , respectively. We assume that 
A = (ci> 2 - ~ (u>2 - w 3 ) - uj c > \g\, |fij|,7e where 

7" 1 is the lifetime of |2). Under these conditions the 
excited state can be adiabatically eliminated, leading to 
two-photon Raman transitions between |1) and |3) with 
the two-photon Rabi frequency x( x ) = gu(x)£li (x) / A 
where u(x) is the mode function of the cavity field. 

We assume that the atoms and the weakly bound 
molecules in state |1) are confined by a trapping po- 
tential inside the optical cavity. The state to be am- 
plified, |3), may also be confined by a trapping potential 
or molecules in that state may pass through the cavity as 
a wave packet ^3 • Here we only treat explicitly the case 
of molecules in |3) being in a stationary state of their 
center-of-mass motion. At zero temperature we can use 
a single mode approximation for the atomic and molec- 
ular fields with the center-of-mass wave functions (j> n (x) 
and energies e n for \n = 0, 1, 2, 3). The Hamiltonian, in a 
rotating frame in which the energy of |3) is given by the 
two-photon detuning, 8 = (1^3 + 63) — (u>i +ei) — (wi — oj c ), 
consists of three terms, H = 5b\b 3 + H\ 3 + Hoi, where 

#13 = \x'b\b 3 a + h.c.,H 01 = a'blboboe 1 ^ + h.c. 

Here, -H13 represents the two-photon transitions between 
|1) and |3) with x' — — J d 3 x4>l(x)x(x)(j) 3 (x) while H m 
represents the Feshbach resonance coupling of states |0) 
and |1) with a' = a J cPxtfii (x)*<Pq(x) and lu[ = u>i + 
(ei — 2eo). The operators bi are bosonic annihilation op- 
erators for atoms or molecules in state |z) and a is the 



annihilation operator for photons in the cavity field. 

The cavity field is damped at a rate k with the damping 
being described by the Born-Markov master equation for 
a damped harmonic oscillator. The time evolution of the 
total density operator, w(t), is then given by dw/dt = 
—i[H, w] + k (2awa) — a) aw — wa'aj. In the bad cavity 
limit where k 3> Ix'IVAo, |^|, where Ng is the maximum 
number of molecules in |1) or |3), the cavity field can be 
adiabatically eliminated to give a master e quat ion for the 
reduced density operator, p — Tr cav i ty [u>] , [15j 

p = -i[H 01 +5blb 3 , p}+ K \/3\ 2 (bibtphbl - kblhblp + h.c?) 

(1) 

and (3 = x' /^ K ^ 1- The term proportional to n\f3\ 2 in 
Eq. Q represents the amplification of molecules in |3). 
Only the |1) — > |3) transition occurs because the photons 
emitted into the cavity are lost before they can be reab- 
sorbed. The first term in Eq.Q, [Hoi,p], represents the 
oscillatory conversion of pairs of atoms into molecules. 
For k |/3 1 2 > a' the coherent oscillations between atom 
pairs and quasi-bound molecules are suppressed and the 
conversion of atoms into molecules dominates over the 
reverse process. 

Equation JQ) can be used to derive the equations of 
motion for the expectation values of the amplitudes, 
(bi) = Tr[bip], and the populations, (hi) = Tr[hip] where 
hi = SjSi, 

(A 3 ) = 2 K |/3| 2 ((l+n 3 )ni) 
(Si) = -«|/3| 2 ((l + n 3 )6i) -la'CbSvV^ 1 
(Ai) = -2«|/3| 2 ((l + n 3 )«i) -o l \i{b b b\)e^ t + c.c.) 
(bok) = -2za'((l + 2n )6 1 )e-^ 4 

and (S3) = -i8(b 3 ) + K,\/3\ 2 (hib 3 ). Note that it is the 
pairing field, (bobo), rather than the the atomic field, (60), 
that drives the formation of molecules in |1). The atomic 
population is determined by conservation of total particle 
number, d((n ) + 2{h 1 ) + 2(h 3 )) jdt = 0. 

These equations do not form a closed set since the ex- 
pectation values of the amplitudes and populations cou- 
ple to the various cross correlations between the modes 
such as (hin 3 ) or (SoSoSj). The simplest approximation 
is to assume that the correlations factorize at all times, 
resulting in the closed set of differential equations, 

h 3 = 2(l+n 8 )n 1 (2) 

Si = -[(l + n 3 )+f/2]& 1 -zaP e^ lT (3) 

hi = -2[(l + n 3 )+T/2}n 1 -a(iP b* 1 e l < :jlT + c.c) (4) 

Po = -2iu(\ + 2no)b x eT iaiT (5) 

n = 2u(iPob\e^ lT + c.c.) (6) 

where rii — (hi), bi = (bi), and Po — (bobo), and the 
dot now indicates a derivative with respect to the di- 
mensionless time r = k|/3| 2 £ with a = a' /k\P\ 2 and 
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u>x = oj[/K\f3\ 2 . We have also introduced a phenomeno- 
logical decay rate, T, for the intermediate state |1) to 
describe losses due to inelastic collisions leading to vi- 
brational decay for ui < or disintegration into pairs of 
atoms for uji > 0. Equations ©-© arc independent of 
the detuning S. This follows from the assumption that 
| <5 1 <C k used to derive Eq. (Q), which implies that the 
two-photon detuning is much less than the linewidth of 
the transition. 

The equation of motion for the amplitude of the final 
state, 63, can be immediately integrated to give 

6 3 (r)=exp[/ dr'H<5/«l/3| 2 + «i(^))]^(0). (7) 
Jo 

Consequently, phase coherent amplification of |3) only 
occurs when there is an initial signal to be amplified, 
&3(0) 7^ 0. Otherwise one only has amplification of the 
population as given by Eq. J5Jl. This remains true even 
for the exact dynamics given by For the remainder 
of this paper we let S = 0. 

Consider first the case u)\ = Y = where |1) is ini- 
tially in the vacuum state and 77,3 <C 1. Then for very 
short times such that the atomic condensate can be 
treated as an undepleted source with Po(0) = 77 (0) = 
2N the population in |3) grows according to 773(7-) = 
exp[2a 2 (2iV) 2 r 3 /3] — 1. This indicates that even for very 
short times the amplification cannot be described via a 
simple linear rate equation of the form dn^/dt = G713 
where G = const. For longer times, b\ and n\ both go 
to zero while 713(7-) saturates at a value that is less than 
its maximal value of 713(0) + N, which implies less than 
100% conversion of atoms into molecules (but approach- 
ing 100% for N — ► 00). Fig. 1 shows the dynamics 
of the populations and amplitudes for the three modes. 
The population and amplitude in |3) grow monotonically 
while states |1) and |0) initially oscillate at a frequency 
- ay/N. 

Figure 2 shows the coherent amplification of 63 and 
the amplification of the population 713, which includes 
both coherent and incoherent contributions, for Po(0) — 
no(0) = 2N = 500, a = 1, and various values of u)i and 
f. As long as u>\ < a\J~N, the |0) — * |1) transition is 
resonant and a finite u)\ has no significant effect on the 
dynamics. For u>x 3> a^/N the transition is off-resonant 
and this significantly slows the rate at which population 
builds up in |3). However, the value at which the pop- 
ulation in |3) saturates is unaffected by <Jj\. Decay from 
the intermediate state only has a significant effect on the 
population and amplitude in |3) for r/2 > 1 + 713(0). In 
this limit the decay out of |1) dominates over the initial 
stimulated scattering from |1) to |3). The fact that the 
dynamics are independent of S (for \S\ <C n) and only 
weakly dependent on ui\ indicate that mean-field energy 
shifts due to two-body interactions between atoms and 
molecules will have a negligible impact on the amplifica- 
tion of 1 3). 
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FIG. 1: (a) Populations of the three modes: ri3 (solid line), ni 
(dotted line), and no/2 (dashed dot line), (b) Amplitudes of 
the three modes: | Z>3 1 (solid line), (dotted line), y // \Po\/2 
(dashed dot line). In both plots, a — 10, wi = T = 0, &i(0) = 
m(0) = 0, n 3 (0) = 63 (0) 2 = 10, and n (0) = P (0) = 2N = 
200. 
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FIG. 2: (a) Population, 713(7"), and (b) amplitude, 63 (t), 
of the ground state molecular field for the initial condi- 
tions 713 (0) = &s(0) 2 = 10, n (0) = P (0) = 500, and 
m(0) = 6i(0) = 0. Solid line: Qi/a = T = 0; Circles: 
Qi/a = 10 and T = 0; Dashed line: Qi/a = 100 and T = 0; 
Dashed dot line: cJi/a = T — 10; Dotted line: u>i/a = 10 and 
f = 100. 
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FIG. 3: (a) Population of the ground state molecular field, 
713 (t), calculated using Eq. (2) (dashed dot line) and directly 
from the master equation, Eq. (1), (solid line), (b) amplitude 
of ground state molecules, 63(7"), from Eq. (7) (dashed dot 
line) and master equation (solid line) . In both figures we have 
taken the initial amplitudes for state |3) to be Co = ci = 1/V2, 
a = 1, u>i = r = 0, and 2N = 16 to be the initial population 
in |0) with |1) initially in the vacuum state. 



We can compare the solution obtained by integrating 
Eqs. ©-© with the expectation values obtained from 
direct integration of the master equation, for small 
numbers of molecules (~ 10). This is shown in Fig. 3. 
For Eq. JQ) we choose the initial conditions to be a Fock 
state with occupation number 2N for mode |0), and |3) 
to be in an arbitrary superposition of or 1 molecules, 
Co \vacuum) + c\b\ \vacuum) . For the factorized equations 
motion we then have the corresponding initial conditions 
n 3 (0) = | Cl | 2 , 6 3 (0) = c* c x , and n (0) = P o (0) = 2N. 

For short times, when the population and field ampli- 
tude in 1 3) are growing exponentially, the two solutions 



show good agreement. For longer times, r > 1/2, the 
two solutions for 713 and 63 saturate at different values. 
Numerical integration of Eq. Q indicates that the con- 
version of atoms into ground state molecules asymptoti- 
cally approaches 100% asr^oo for all N in contrast to 
the results obtained from Eqs. This difference 

can be explained by noting that 713 (t) and &3(t) saturate 
when no more new population forms in |1). From the 
equation of motion for (hi) one can see that the pop- 
ulation in |1) is driven by the cross-correlation (b\bobo). 
When this goes to zero, no more population is transferred 
from |0) to |1) and as a result n^r) saturates. Figure 4 
shows (Sj&o^o) calculated directly from the master equa- 
tion and b\P$ calculated using Eqs. ©-©■ As one can 
see, b\P$ decays to zero much faster than (&j&o^o)- Con- 



FIG. 4: (feJSoSo) calculated from the master equation (solid 
line) and b\ Pq calculated from factorized equations of motion 
(dashed dot line). All parameters are the same as Fig. 3. 



sequently the buildup of population in 1 1) as given by Eq. 
(0J stops well before that given by the exact quantum dy- 
namics. We note that other cross-correlations involving 
the intermediate state such as (n3?ii) exhibit similar be- 
havior with the uncorrelated factorized values obtained 
from Eqs. ©-© decaying to zero much faster than the 
exact values obtained from the master equation. 

In conclusion we have analyzed a model for a coherent 
matter-wave amplifier of ground-state molecules in an 
optical cavity by analyzing semiclassical rate equations 
and comparing them to the exact quantum dynamics. 

This work is supported in part by the US Office of 
Naval Research, the NSF, the US Army Research Office, 
NASA, and the Joint Services Optics Program. 
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